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Abstract. In spite of the great progress we have seen in recent years in the derivation 
of nuclear forces from chiral effective field theory (EFT), some important issues are still 
unresolved. In this contribution, wo discuss the open problems which have particular 
relevance for microscopic nuclear structure, namely, the proper renormalization of 
chiral nuclear potentials and sub-leading many-body forces. 
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1. Introduction 

The fundamental goal of nuclear structure physics is to understand the properties of 
atomic nuclei and their reactions in terms of the basic forces between the constituents. 
During the past half century, a large variety of phenomenological forces has been 
developed and applied in microscopic nuclear structure calculations with some success. 
But in the long run phenomenology is not good enough and, ultimately, we need nuclear 
interactions that are based upon proper theory. Since the nuclear force is a manifestation 
of strong interactions, any serious derivation has to start from quantum chromodynamics 
(QCD). However, the well-known problem with QCD is that it is non-perturbative in 
the low-energy regime characteristic for nuclear physics. For many years this fact was 
perceived as the great obstacle for a derivation of nuclear forces from QCD — impossible 
to overcome except by lattice QCD. The effective field theory (EFT) concept has shown 
the way out of this dilemma. One has to realize that the scenario of low-energy QCD is 
characterized by pions and nucleons interacting via a force governed by spontaneously 
broken approximate chiral symmetry. This chiral EFT allows for a systematic low- 
momentum expansion known as chiral perturbation theory (ChPT) [T]. Contributions 
are analyzed in terms of powers of small external momenta over the large scale: {Q/A^y, 
where Q is generic for an external momentum (nucleon three-momentum or pion four- 
momentum) or pion mass and ^ 1 GeV is the chiral symmetry breaking scale ('hard 
scale'). The early applications of ChPT focused on systems like tttt [2] and nN [3], where 
the Goldstone-boson character of the pion guarantees that the expansion converges. 

The past 15 years have also seen great progress in applying ChPT to nuclear forces 
[lElElEllHlElIinillllinillSlElttSlESlinillH]. As a result, nucleon-nucleon (A^A^) 
potentials of high precision have been constructed, which are based upon ChPT carried 
to next-to-next-to-next-to-leading order (N'^LO) [151 [HI [H], and applied in nuclear 
structure calculations with great success. 

However, in spite of this progress, we are not done. Due to the complexity of the 
nuclear force issue, there are still many subtle and not so subtle open problems. We will 
not list and discuss all of them, but instead just focus on what is relevant to nuclear 
structure physics. In this regard, there are two important issues that need our attention: 

• The proper renormalization of chiral nuclear potentials and 

• Subleading chiral few-nucleon forces. 

To set the stage, we will give first, in section 2, a brief overview of chiral nuclear forces 
and ChPT. The main two open issues are discussed in section 3 and 4. The article is 
concluded in section 5. 

2. Nuclear forces, chiral perturbation theory, and power counting 

Effective field theories (EFTs) are defined in terms of effective Langrangians which 
are given by an infinite series of terms with increasing number of derivatives and/or 
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nucleoli fields, with the dependence of each term on the pion field prescribed by the 
rules of broken chiral symmetry. Applying this Lagrangian to a particular process, an 
unlimited number of Feynman graphs can be generated. Therefore, we need a scheme 
that makes the theory manageable and calculabel. This scheme which tells us how to 
distinguish between large (important) and small (unimportant) contributions is chiral 
perturbation theory (ChPT), and determining the power i> of the expansion has become 
known as power counting. 

Nuclear potentials are defined as sets of irreducible graphs up to a given order. 
The power z/ of a few-nucleon diagram involving A nucleons is given in terms of naive 
dimensional analysis by: 

z/ = -2 + 2A-2C + 2L + ^Ai, (1) 

i 

with 

A.^rf. + ^-2, (2) 

where C denotes the number of separately connected pieces and L the number of loops in 
the diagram; di is the number of derivatives or pion-mass insertions and the number of 
nucleon fields (nucleon legs) involved in vertex i; the sum runs over all vertices contained 
in the diagram under consideration. Note that Aj > for all interactions allowed by 
chiral symmetry. For an irreducible NN diagram ( "two-nucleon force", A = 2, C = 1), 
([T]) collapses to 

u = 2L + Y.Ai. (3) 

i 

The power formula ([T]) allows to predict the leading orders of multi-nucleon forces. 
Consider a m-nucleon irreducibly connected diagram (m-nucleon force) in an A-nucleon 
system [m < A). The number of separately connected pieces isC = A — m + 1. Inserting 
this into ([T| together with L = and J2i Aj = yields u = 2m — 4. Thus, two-nucleon 
forces (m = 2) start at z/ = 0, three-nucleon forces (m = 3) at z/ = 2 (but they happen 
to cancel at that order), and four-nucleon forces at z/ = 4 (they don't cancel). Thus, 
ChPT provides a straightforward explanation for the empirically known fact that 2NF 
> 3NF > 4NF .... 

In summary, the chief point of the ChPT expansion is that, at a given order u, there 
exists only a finite number of graphs. This is what makes the theory calculable. The 
expression {Q/A^Y'^^ provides a rough estimate of the relative size of the contributions 
left out and, thus, of the accuracy at order z/. In this sense, the theory can be calculated 
to any desired accuracy and has predictive power. 

Chiral perturbation theory and power counting imply that nuclear forces emerge 
as a hierarchy controlled by the power z/ (see [16j for a pedagogial introduction). 

Since 2003, a very quantitative chiral A^A^ potential (at N^LO, u = 4) [TS] exists 
which has been applied successfully in many nuclear structure calculations [191 EOl EH 
[22| [23] . However, questions have been raised concerning the proper renormalization 
of chiral A^A^ potentials, which is why we will look into this issue in the next section. 
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Moreover, there are also some open problems in the few-nucleon-force sector, which will 
be discussed in section HI 

3. Renormalization of chiral nuclear forces 

3.1. The chiral NN potential 

In terms of naive dimensional analysis or "Weinberg counting" , the various orders of 
the irreducible graphs which define the chiral NN potential are given by: 

V,o =Vr + Vrj (4) 
^NLo =V^o + vi'^ + vl!^ + vS^ (5) 

^NNLO = ^NLO + + Vff (6) 

V^N3LO = VS,NLO + ^ + Vl:^ + + V^£^ (7) 

where the superscript denotes the order u of the low-momentum expansion. LO stands 
for leading order, NLO for next-to-leading order, etc.. Contact potentials carry the 
subscript "ct" and pion-exchange potentials can be identified by an obvious subscript. 
The one-pion exchange (IPE) potential reads 

ViniP ,P) = - jj^ Ti ■ T2 ^ , (8 

4/2 g2 + ^2 

where p ' and p designate the final and initial nucleon momenta in the center-of-mass 
system and q = p'—pis the momentum transfer; ai 2 and Ti 2 are the spin and isospin 
operators of nucleon 1 and 2; g^, fn, and denote axial-vector coupling constant, the 
pion decay constant, and the pion mass, respectively. Since higher order corrections 
contribute only to mass and coupling constant renormalizations and since, on shell, 
there are no relativistic corrections, the on-shell IPE has the form ([s]) up to all orders. 

Multi-pion exchange, which starts at NLO and continues through all higher orders, 
involves divergent loop integrals that need to be regularized. An elegant way to 
do this is dimensional regularization which (besides the main nonpolynomial result) 
typically generates polynomial terms with coefficients that are, in part, infinite or 
scale dependent p. One purpose of the contacts is to absorb all infinities and scale 
dependencies and make sure that the final result is finite and scale independent. 
This is the renormalization of the perturbatively calculated A^A^ amplitude (which, 
by definition, is the 'WA^ potential"). It is very similar to what is done in the ChPT 
calculations of tttt and vrA^ scattering, namely, a renormalization order by order, which is 
the method of choice for any EFT. Thus, up to this point, the calculation fully meets the 
standards of an EFT and there are no problems. The perturbative A^A^ amplitude can 
be used to make model independent predictions for peripheral partial waves [9l HOj [H] . 

3.2. Nonperturbative applications of the NN potential 



For calculations of the structure of nuclear few and many-body systems, the lower 
partial waves are the most important ones. The fact that in 5* waves we have large 
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scattering lengths and shallow (quasi) bound states indicates that these waves need to 
be treated nonperturbatively. Following Weinberg's prescription [1] , this is accomplished 
by inserting the potential V into the Lippmann-Schwinger (LS) equation: 

T{p',p) = V{p',p) + f £p" V{p',p") 2 _^,^2 I ■ T{p",p) , (9) 

where denotes the nucleon mass. 

In general, the integral in the LS equation is divergent and needs to be regularized. 
One way to do this is by multiplying V with a regulator function 

Vip \p)^ V{p p) e-(P'/^)'" e-(*'/'^)'" . (10) 

Typical choices for the cutoff parameter A that appears in the regulator are A ^ 
0.5 GeV < A^ ^ 1 GeV. 

It is pretty obvious that results for the T-matrix may depend sensitively on the 
regulator and its cutoff parameter. This is acceptable if one wishes to build models. 
For example, the meson models of the past [Ml ES] always depended sensitively on the 
choices for the cutoff parameters which, in fact, were important for the fit of the NN 
data. However, the EFT approach wishes to be fundamental in nature and not just 
another model. 

In field theories, divergent integrals are not uncommon and methods have been 
developed for how to deal with them. One regulates the integrals and then removes the 
dependence on the regularization parameters (scales, cutoffs) by renormalization. In the 
end, the theory and its predictions do not depend on cutoffs or renormalization scales. 

So-called renormalizable quantum field theories, like QED, have essentially one set 
of prescriptions that takes care of renormalization through all orders. In contrast, EFTs 
are renormalized order by order. 

As discussed, the renormalization of perturbative EFT calculations is not a problem. 
The problem is nonperturbative renormalization. This problem typically occurs in 
nuclear EFT because nuclear physics is characterized by bound states which are 
nonperturbative in nature. EFT power counting may be different for nonperturbative 
processes as compared to perturbative ones. Such difference may be caused by the 
infrared enhancement of the reducible diagrams generated in the LS equation. 

Weinberg's implicit assumption [H |26] was that the counterterms introduced to 
renormalize the perturbatively calculated potential, based upon naive dimensional 
analysis ("Weinberg counting"), are also sufficient to renormalize the nonperturbative 
resummation of the potential in the LS equation. In 1996, Kaplan, Savage, and 
Wise (KSW) [27] pointed out that there are problems with the Weinberg scheme if 
the LS equation is renormalized by minimally-subtracted dimensional regularization. 
This criticism resulted in a flurry of publications on the renormalization of the 
nonperturbative A^A^ problem [28l[29llMlEll3a[33l[Sll35l|36l[371[^ 
The literature is too comprehensive to discuss all contributions. Let us just mention 
some of the work that has particular relevance for our present discussion. 
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If the potential V consists of contact terms only (a.k.a. pion-less theory), then the 
nonperturbative summation ^ can be performed analytically and the power counting 
is explicit. However, when pion exchange is included, then ^ can be solved only 
numerically and the power counting is less transparent. Perturbative ladder diagrams 
of arbitrarily high order, where the rungs of the ladder represent a potential made 
up from irreducible pion exchange, suggest that an infinite number of counterterms is 
needed to achieve cutoff independence for all the terms of increasing order generated by 
the iterations. For that reason, Kaplan, Savage, and Wise (KSW) [27] proposed to sum 
the leading-order contact interaction to all orders (analytically) and to add higher-order 
contacts and pion exchange perturbatively up to the given order. Unfortunately, it 
turned out that the order by order convergence of IPE is poor in the ^^i-^Z^i state [28]. 
The failure was triggered by the singularity of the IPE tensor force when iterated 
to second order. Therefore, KSW counting is no longer taken into consideration (see, 
however, [12]). A balanced discussion of possible solutions can be found in [32] . 

Some researchers decided to take a second look at Weinberg's original proposal. A 
systematic investigation of Weinberg counting in leading order has been conducted by 
Nogga, Timmermans, and van Kolck [31| in momentum space, and by Valderrama and 
Arriola at LO and higher orders in configuration space [331 ESI [36] . A comprehensive 
discussion of both approaches and their equivalence can be found in [311 HSj. 

The LO A^A^ potential is given in Q and consists of IPE plus two nonderivative 
contact terms that contribute only in 5* waves. Nogga et al find that the given 
counterterms renormalize the S waves (i.e., stable results are obtained for A — )• oo) 
and the naively expected infinite number of counterterms is not needed. This means 
that Weinberg power counting does actually work in S waves at LO (ignoring the 
dependence of the contact interaction discussed in Refs. [271 [32]). However, there are 
problems with a particular class of higher partial waves, namely those in which the tensor 
force from IPE is attractive. The first few cases of this kind of low angular momentum 
are ^Pq, ^-^2, and ^-D2, which need a counterterm for cutoff independence. The leading 
order (nonderivative) counterterms do not contribute in P and higher waves, which is 
why Weinberg counting fails in these cases. But the second order contact potential 
provides counterterms for P waves. Therefore, the promotion of, particularly, the ^Pq 
and ^P2 contacts from NLO to LO would fix the problem in P waves. To take care of the 

problem, a N^LO contact, i.e. a term from V^t^"*, needs to be promoted to LO. Partial 
waves with orbital angular momentum L > 3 may be calculated in Born approximation 
with sufficient accuracy and, therefore, do not pose renormalization problems. In this 
way, one arrives at a scheme of 'modified Weinberg counting' [31] for the leading order 
two-nucleon interaction. 

3.3. Renormalization beyond leading order 

As discussed, for a quantitative chiral A^A^ potential one needs to advance all the way to 
N^LO. Thus, the renormalization issue needs to be discussed beyond LO. Naively, the 
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most perfect renormalization procedure is the one where the cutoff parameter A is carried 
to infinity while stable results are maintained. This was done successfully at LO in the 
work by Nogga et al [3l] described above. At NNLO, the infinite-cutoff renormalization 
procedure has been investigated in [ID] for partial waves with total angular momentum 
J < 1 and in for all partial waves with J < 5. At N^LO, an investigation of the 
^5*0 state exists [39]. From all of these works, it is evident that no counter term is 
effective in partial-waves with short-range repulsion and only a single counter term can 
effectively be used in partial-waves with short-range attraction. Thus, for the A — > oo 
renormalization prescription, even at N^LO, there is either one or no counter term per 
partial- wave state. This is inconsistent with any reasonable power-counting scheme and, 
therefore, defies the principals of an EFT. 

A possible way out of this dilemma was proposed already in [31] and reiterated in 
a recent paper by Long and van Kolck [H]. In the latter reference, the authors examine 
the renormalization of an attractive potential perturbed by a correction. 
Generalizing their findings, they come to the conclusion that, for any attractive 1/r"' 
potential (with n > 2), partial waves with low angular momentum L must be summed to 
all orders and one contact term is needed for each L to renormalize the LO contribution. 
However, there exists an angular momentum Lp {Lp ^ 3 for the nuclear case, cf. 
[31]), above which the leading order can be calculated perturbatively. In short, naive 
dimensional analysis (NDA) does not apply at LO below Lp. However, once this failure 
of NDA is corrected at LO, higher order corrections can be added in perturbation theory 
using counterterm that follow NDA [H]. 

Reference [H] used just a toy model and, therefore, a full investigation using the 
chiral expansion is needed to answer the question if this renormalization approach will 
work for the realistic nuclear force. A first calculation of this kind for the S waves was 
recently performed by Valderrama [H]. The author renormalizes the LO interaction 
nonperturbatively with A — oo and then uses the LO distorted wave to calculate the 
2PE contributions at NLO and NNLO perturbatively. It turns out that perturbative 
renormalizability requires the introduction of three counterterms in ^5*0 and six in the 
coupled ^Si Di channels. Thus, the number of counterterms required in this scheme 
is larger than in the Weinberg scheme, which reduces the predictive power. For a final 
evaluation of this approach, also the results for P and D waves are needed, which are 
not yet available. 

However, even if such a project turns out to be successful for A^A^ scattering, there 
is doubt if the interaction generated in this approach is of any use for applications in 
nuclear few- and many-body problems. In applications, one would first have to solve 
the many-body problem with the renormalized LO interaction, and then add higher 
order corrections in perturbation theory. However, it was shown in a recent paper [13] 
that the renormalized LO interaction is characterized by a very large tensor force from 
IPE. This is no surprise since LO is renormalized with A — )■ cxd implying that the IPE, 
particulary its tensor force, is totally uncut. As a consequence of this, the wound integral 
in nuclear matter, k, comes out to be about 40%. The hole-line and coupled cluster 
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expansions are know to converge ~ with n the number of hole-hnes or particles 
per cluster |l6l H?! HH] . For conventional nuclear forces, the wound integral is typically 
between 5 and 10% and the inclusion of three-body clusters (or three hole-lines) are 
needed to obtain converged results in the many-body system [201 ISSl HI] . Thus, if the 
wound integral is 40%, probably, up to six hole-lines need to be included for approximate 
convergence. Such calculations are not feasible even with the most powerful computers 
of today and will not be feasible any time soon. Therefore, even if the renormalization 
procedure proposed in [H] will work for A^A^ scattering, the interaction produced will 
be highly impractical (to say the least) in applications in few- and many-body problems 
because of convergence problems with the many-body energy and wave functions. 

3.4- Back to the beginnings 

The various problems with the renormalization procedures discussed above may have 
a simple common reason: An EFT that has validity only for momenta Q < is 
applied such that momenta Q ^ A^ are heavily involved (because the regulator cutoff 
A — )■ oo). A recent paper by Epelbaum and Gegelia |SU] illustrates the point: The 
authors construct an exactly solvable toy-model that simulates a pionful EFT and yields 
finite results for A — )■ cxd. However, as it turns out, these finite results are incompatible 
with the underlying EFT, while for cutoffs in the order of the hard scale consistency is 
maintained. In simple terms, the point to realize is this: // an EFT calculation produces 
(accidentally) a finite result for A — )■ oo, then that does not automatically imply that 
this result is also right. 

This matter is further elucidated in the lectures by Lepage of 1997 ^T\. Lepage 
points out that it makes little sense to take the momentum cutoff beyond the range of 
validity of the effective theory. By assumption, our data involves energies that are too 
low — wave lengths that are too long — to probe the true structure of the theory at very 
short distances. When one goes beyond the hard-scale of the theory, structures are seen 
that are almost certainly wrong. Thus, results cannot improve and, in fact, they may 
degrade or, in more extreme cases, the theory may become unstable or untunable. In 
fact, in the A^A^ case, this is what is happening in several partial waves (as reported 
above). Therefore, Lepage suggests to take the following three steps when building an 
effective theory: 

(i) Incorporate the correct long-range behavior: The long-range behavior of the 
underlying theory must be known, and it must be built into the effective theory. In 
the case of nuclear forces, the long-range theory is, of course, well known and given 
by one- and multi-pion exchanges. 

(ii) Introduce an ultraviolet cutoff to exclude high-momentum states, or, equivalent, to 
soften the short-distance behavior: The cutoff has two effects: First it excludes high- 
momentum states, which are sensitive to the unknown short-distance dynamics; 
only states that we understand are retained. Second it makes all interactions regular 
at r = 0, thereby avoiding the infinities that beset the naive approach. 
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(iii) Add local correction terms (also known as contact or counter terms) to the effective 
Hamiltonian. These mimic the effects of the high-momentum states excluded by 
the cutoff introduced in the previous step. In the meson-exchange picture, the 
short-range nuclear force is described by heavy meson exchange, like the p(770) 
and uj{782). However, at low energy, such structures are not resolved. Since we 
must include contact terms anyhow, it is most efficient to use them to account for 
any heavy-meson exchange as well. The correction terms systematically remove 
dependence on the cutoff. 

A first investigation in the above spirit has been conducted by Epelbaum and 
MeiBner in 2006. The authors stress that there is no point in taking the cutoff 
A beyond the breakdown scale of the EFT, A-^ ^ mp ^ 1 GeV, since the error of 
the calculation is not expected to decrease in that regime. Any value for the cutoff 
parameter A is acceptable if the error associated with its finite value is within the 
theoretical uncertainty at the given order. The authors conduct an investigation at LO 
(including only the counter terms implied by Weinberg counting) and find that, starting 
from A 3 fm~^, the error in the A^A^ phase shifts due to keeping A finite stays within 
the theoretical uncertainty at LO. 

3. 5. Bringing the renormalization business to a finish 

Crucial for an EFT are regulator independence (within the range of validity of 
the EFT) and a power counting scheme that allows for order-by-order improvement 
with decreasing truncation error. The purpose of renormalization is to achieve this 
regulator independence while maintaining a functional power counting scheme. After the 
comprehensive tries and errors of the past, it appears that there are two renormalization 
schemes which have the potential to achieve the above goals and, therefore, should be 
investigated systematically in the near future — with the hope to bring the tiresome 
renormalization issue finally to conclusion. 

In scheme one, the LO calculation is conducted nonperturbatively (with A — )■ oo 
as in [31]) and subleading orders are added perturbatively in distorted wave Born 
approximation. As mentioned above, Valderrama has started this in 5" waves [H], but 
results in higher partial waves are needed to fully assess this approach. Even though 
at this early stage any judgement is speculative, we take the liberty to predict that 
this approach will be only of limited success and utility — for the following reasons. 
First, it will probably require about twice as many counterterms as Weinberg counting 
and, therefore, will have less predictive power. Second, this scheme may converge 
badly, because the largest portion of the nuclear force, namely, the intermediate-range 
attraction appears at NNLO. Third, as discussed in [15], this force may be problematic 
(and, therefore, impractical) in applications in nuclear few- and many-body systems, 
because of a pathologically strong tensor force that will cause bad convergence of energy 
and wave functions. Finally, in the work that has been conducted so far within this 
scheme by Valderrama, it is found that only rather soft cutoffs can be used. 
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Figure 1. The three-nucleon force at NNLO. From left to right: 2PE, IPE, and 
contact diagrams. Sohd hnes represent nucleons and dashed Unes pions. Smah dots 
denote vertices with = and large solid dots are = 1. 

The latter point (namely, soft cutoffs) suggests that one may then as well conduct 
the calculation nonperturbatively at all orders (up to N'^LO) using Weinberg counting, 
which is no problem with soft cutoffs. This is scheme two that we propose to investigate 
systematically. In the spirit of Lepage, the cutoff independence should be examined 
for cutoffs below the hard scale and not beyond. Ranges of cutoff independence within 
the theoretical error are to be identified using 'Lepage plots' [5l]. A very systematic 
investigation of this kind does not exist at this time and is therefore needed, once and for 
all and for principal reasons. Comprehensive circumstantial evidence from the numerous 
chiral A^A^ potentials constructed over the past decade [HI [TT], [131 UHl [161 IE! may be 
perceived as an indication that this investigation will most likely be a success. If so, 
then the results of such investigation will resolve the renormalization issue and remove 
concerns about existing chiral A^A^ potentials [151 [IZj, which are currently applied in 
microscopic nuclear structure physics with great success. 

4. Few-nucleon forces and what is missing 

We will now discuss the other issue we perceive as unfinished and important, namely, 
subleading chiral few-nucleon forces. 

Nuclear three-body forces in ChPT were initially discussed by Weinberg [5]. The 
3NF at NNLO, was derived by van Kolck [7| and applied, for the first time, in nucleon- 
deuteron scattering by Epelbaum et al [22] • The leading 4NF (at N^LO) was recently 
constructed by Epelbaum [53j and found to contribute in the order of 0.1 MeV to the "^He 
binding energy (total ^He binding energy: 28.3 MeV) in a preliminary calculation [M] . 
confirming the traditional assumption that 4NF are essentially negligible. Therefore, 
the focus is on 3NF. 

For a 3NF, we have A = ?> and C = 1 and, thus, ([T| implies for 3NF 

V = 2 + 2L + Y,/\,. (11) 

i 

We will use this equation to analyze 3NF contributions order by order. The lowest 
possible power is obviously v = 2 (NLO), which is obtained for no loops (L = 0) and 
only leading vertices (I]j Aj = 0). This 3NF happens to vanish [5l |55l |56]. The first 
non-vanishing 3NF occurs at NNLO. 
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Figure 2. The 3NF at N^LO: Two examples of one-loop graphs. Notation as in 
figure [l] 



4.1. The 3NF at NNLO 

The power z/ = 3 (NNLO) is obtained when there are no loops (L = 0) and Y^i^i = 1, 
i.e., Aj = 1 for one vertex while Aj = for all other vertices. There are three topologies 
which fulfill this condition, known as the two-pion exchange (2PE), IPE, and contact 
graphs (figure [T| . 

The 3NF at NNLO has been derived (without the I/Ma? corrections) [3 [52] and 
applied in calculations of few-nucleon reactions EZl EHl EH] , structure of light- and 
medium-mass nuclei |2Tl [601 EI] , and nuclear and neutron matter [23l [62] with a fair 
deal of success. However, the famous 'Ay puzzle' of nucleon-deuteron scattering is 
not solved [521 ISZ]) the even bigger problem with the analyzing power in p-^He 
scattering [631 El] will certainly not be fixed at this order. Furthermore, the spectra of 
light nuclei leave room for improvement |21j . 

We note that there are further 3NF contributions at NNLO, namely, the 1/Mn 
corrections of the NLO 3NF diagrams. Part of these corrections have been calculated 
by Coon and Friar in 1986 [56]. These contributions are believed to be very small. 

In summary, because of various unresolved problems in microscopic nuclear 
structure, the 3NF beyond NNLO is very much in need. In fact, it is no exaggeration 
to state that the 3NF at sub-leading orders is presently one of the most important 
outstanding issues in the chiral EFT approach to nuclear forces. 



4.2. The 3NF at I^LO 



According to (11), the value z/ = 4, which corresponds to N LO, is obtained for the 
following classes of diagrams. 



3NF loop diagrams at ISPLO. For this group of graphs, we have L = 1 and, therefore, 
all Aj have to be zero to ensure u = 4. Thus, these one-loop 3NF diagrams can 
include only leading order vertices, the parameters of which are fixed from ttN and A^A^ 
analysis. We show two samples of this very large class of diagrams in figure [2] One sub- 
group of these diagrams ( "27r exchange graphs" ) has been calculated by Ishikawa and 
Robilotta [65], and two other topologies (27r-l7r and ring diagrams) have been evaluated 
by the Bonn-Jiilich group [66]. The remaining topologies, which involve a leading order 
four-nucleon contact term (e.g., second diagram of figure |2|), are under construction by 
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the Bonn-Jiilich group. The N^LO 27r-exchange 3NF has been apphed in the calculation 
of nucleon-deuteron observables in [65] producing very small effects. 

The smallness of the 27r loop 3NF at N^LO is not unexpected. It is consistent with 
experience with corresponding 2NF diagrams: the NLO 2PE contribution to the NN 
potential, which involves one loop and only leading vertices, is also relatively small. 

By the same token, one may expect that also all the other N^LO 3NF loop topologies 
will produce only small effects. 

3NF tree diagrams at ISPLO. The order i/ = 4 is also obtained for the combination 
L = (no loops) and X^i^j = 2. Thus, either two vertices have to carry Aj = 1 or one 
vertex has to be of the Aj = 2 kind, while all other vertices are Aj = 0. This is achieved if 
in the NNLO 3NF graphs of figure [T] the power of one vertex is raised by one. The latter 
happens if a relativistic l/M^r correction is applied. A closer inspection reveals that 
all l/Mjsf corrections of the NNLO 3NF vanish and the first non- vanishing corrections 
are proportional to and appear at N^LO. However, there are non-vanishing 

corrections of the NLO 3NF and there are so-called drift corrections [HZ] which contribute 
at N^LO (some drift corrections are claimed to contribute even at NLO [67]). We do 
not expect these contributions to be sizable. Moreover, there are contributions from 
the Aj = 2 Lagrangian [6H] proportional to the low-energy constants di. As it turns 
out, these terms have at least one time-derivative, which causes them to be Q/Mn 
suppressed and demoted to N^LO. 

Thus, besides some minor corrections, there are no tree contributions to the 

3NF at N^LO. 

Summarizing the entire N^LO 3NF contribution: For the reasons discussed, we 
anticipate that this 3NF is weak and will not solve any of the outstanding problems. In 
view of this expectation, we have to look for more sizable 3NF contributions elsewhere. 

4.3. The 3NF at N^LO of the A-less theory 

The obvious step to take is to proceed to the next order, N'^LO or z/ = 5, of the A-less 
theory which is the one we have silently assumed so far. (The A-full theory will be 
introduced and discussed below.) Some of the tree diagrams that appear at this order 
were mentioned already: the corrections of the NNLO 3NF and the trees with 

one di vertex which are l/M^ suppressed. Because of the suppression factors, we do 
not expect sizable effects from these graphs. Moreover, there are also tree diagrams 
with one vertex from the Aj = 3 vrA^ Lagrangian [^ [70] proportional to the LECs Cj. 
Because of the high dimension of these vertices and assuming reasonable convergence, 
we do not anticipate much from these trees either. 

However, we believe that the loop contributions that occur at this order are truly 
important. They are obtained by replacing in the N^LO loops (figure [2]) one vertex by 
a Aj = 1 vertex [with LEG Cj]. We show one symbolic example of this large group of 
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a 



Figure 3. (a) One-loop 3NF at N^LO of the A-less theory, (b) Corresponding diagram 
of the A-fuh theory which contributes at N'^LO. Double lines represent A isobars; other 
notation as in figure [T] 



diagrams in figure [3|^a). This 3NF is presumably large and, thus, what we are looking 
for. 

The reasons, why these graphs are large, can be argued as follows. Corresponding 
2NF diagrams are the three-pion exchange (3PE) contributions to the A^A^ interaction. 
In analogy to Figs. |2] and |3](a), there are 3PE 2NF diagrams with only leading vertices 
and the ones with one (sub- leading) Cj vertex (and the rest leading). These diagrams 
have been evaluated by Kaiser in Refs. [71] and [72], respectively. Kaiser finds that 
the 3PE contributions with one sub-leading vertex are about an order magnitude larger 
then the leading ones. 

4.4- N^LO 3NF contributions in the A-full theory 

The above considerations indicate that the A-less theory exhibits, in some cases, a bad 
convergence pattern. The reason for the unnaturally strong subleading contributions 
are the large values of the Aj = 1 LECs, q. The large values can be explained in terms 
of resonance saturation [73] . The A(1232)-resonance contributes considerably to C3 and 
C4. The explicit inclusion of the A takes strength out of these LECs and moves this 
strength to a lower order, thus improving the convergence [TH O [TOl [751 ES]. Figure |3] 
illustrates this fact for the 3NF under consideration: the diagram of the A-less theory 
shown in (a) is (largely) equivalent to diagram (b) which includes one A excitation. 
Note, however, that diagram (a) is N^LO, while diagram (b) is N^LO. Moreover, there 
are further N^LO one-loop diagrams with two and three A excitations, which correspond 
to diagrams of order N^LO and N^LO, respectively, in the A-less theory. 

This consideration clearly shows that the inclusion of A degrees of freedom in 
chiral EFT makes the calculation of sizable higher-order 3NF contributions much more 
efficient. 

4-5. Summarizing the open 3NF business 

To make a complicated story short, this is the bottom line concerning 3NF: 

• The chiral 3NF at NNLO is insufficient. Additional sizable 3NF contributions are 
needed. 
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• The chiral 3NF at N^LO (in the A-less theory) most hkely does not produce sizable 
contributions. 

• Sizable contributions are expected from one-loop 3NF diagrams at N^LO of the 
A-less or N^LO of the A-fuU theory (figure [s]). These 3NF contributions may turn 
out to be the missing pieces in the 3NF puzzle and have the potential to solve the 
outstanding problems in microscopic nuclear structure^ 

5. Conclusions and Outlook 

The past 15 years have seen great progress in our understanding of nuclear forces in terms 
of low-energy QCD. Key to this development was the realization that low-energy QCD 
is equivalent to an effective field theory (EFT) which allows for a perturbative expansion 
that has become know as chiral perturbation theory (ChPT). In this framework, two- 
and many-body forces emerge on an equal footing and the empirical fact that nuclear 
many-body forces are substantially weaker than the two-nucleon force is explained 
automatically. 

In spite of the great progress and success of the past 15 years, there are still 
some unresolved issues that will need our attention in the near future. One problem 
is the proper renormalization of the chiral two- and many-nucleon potentials. We have 
discussed this issue in section 3, where we also spelled out the systematic work that 
needs to be done to resolve the problems. 

The other unfinished business are the few-nucleon forces beyond NNLO ("sub- 
leading few-nucleon forces") which are needed to hopefully resolve some important 
outstanding nuclear structure problems. We believe that we identified correctly where 
these forces will emerge within the systematic scheme of ChPT. 

If the open issues discussed in this paper will be resolved within the next few years, 
then, after 80 years of desperate struggle, we may finally claim that the nuclear force 
problem is essentially solved. The greatest beneficiary of such progress will be the field 
of ab initio nuclear structure physics. 
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